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Abstract 

We provide explicit formulas for the Green function of an elliptic PDE in the infinite strip 
and the half-plane. They are expressed in elementary and special functions. Proofs of uniqueness 
and existence are also given. 

1 Motivation 

This research is motivated by the multi-asset American option pricing problem(see e.g. [2]). Us¬ 
ing the standard arbitrage theory framework it can be shown that the option prices are bounded 
by solutions of elliptic PDE problems. These problems have specific features: the domains are 
unbounded and the boundary functions are not smooth and do not vanish at infinity. We reduce 
the original PDE determining the price to a static Klein-Gordon equation (SKGE) in an infinite 
strip (for the alternative dual options) and in the half plane (the exchange and basket options). 
The boundary conditions of the corresponding problems are bounded Holder functions. This 
leads to the following questions: 

i. Can we explicitly solve the boundary problems for the SKGE? 

ii. If this is the case, can we construct computer- friendly representations? 

iii. Is the obtained solution classical? 

iv. Is it unique? 

To our knowledge, these aspects of the mentioned problems have not been studied in the liter¬ 
ature. Even the answers to the two last questions seem not be obvious. The majority of the 
known results deal with bounded domains and one can not directly apply available theorems to 
the equations in the infinite strip and in the half-plane. Also, if the boundary conditions are 
not twice differentiable, then the smoothness of the solution is not clear. In this paper we give 
answers to all four questions. Our results are as follows: 

i. The boundary value problem in the infinite strip and in the half-plane allows a closed form 
solution. 

ii. Solutions can be represented in terms of elementary functions or in terms of special func¬ 
tions. 

*1 am grateful to Prof. Yu.M. Kabanov and to Dr. M. Zhitlukhin for important discussions and help 
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iii. For all bounded boundary functions with Holder property the solution is classical. 

iv. The solution is unique. 

We could find in the literature only a few related works. The paper [4] contains closed form 
formulas of boundary problem for the Laplace equation in the infinite strip. In the recent paper 
[B] there is a variety of Green functions for the homogeneous Poisson problems for SKGE in some 
unbounded domains, including the strip and the half-plane. 

The paper is organized as follows. In Section [2] we obtain explicit formulas for the Green 
function in the case of infinite strip and formulate the main theorem on existence, uniqueness 
and closed form of the solution. In Section [3] we address the case of half-plane domain. All proofs 
are given in Section [4j 


2 Problem in the infinite strip 

2.1 Existence and uniqueness theorem 

Let n 1 = lx [0,7r] = {(x, y) € R 2 , y € [0,7r]}. Let LV = AV — r 2 V be an operator acting 
on the twice differentiable functions V = V(x,y ) defined in the interior of TV. Here A is the 
Laplacian and re®. We consider the boundary value problem 


LV(x, y) = 0, (; x , y) £ int n 71 

V(x, 0) = tp(x), 

V(x, 7r) = 0, i£l. 


( 2 . 1 ) 


The following theorem claims that under certain assumptions the problem (2.1) admits a unique 
classical solution and provides an explicit form for it. 

Theorem 2.1. Let H x be the space of Holder func tions of order A > 0 and let be a bounded 
function from H x .Then the solution of the problem (2.1) exists in the classical sense, is unique, 
and allows the representation 


V(x,y)= [ <p(u)G w (x 

JR 


u,y)du , 


( 2 . 2 ) 


with the Green function 


G*(x,y)=5(x)Q(-y) + -J2 

7T 

k =1 


k sin ky p _^ vw ^ P i 
Jk 2 + r 2 


(2.3) 


where S(x) is the Dirac delta-function and 0(x ) = I( o,oo)( x )- 

Corollary 2.1. The Green function G\ for the Laplace equation AH = 0 (i.e. for r = 0) has 
the representation (see W : 


G\[x,y) = 8{x)Q{-y) + 


1 siny 
27 t cosh x — cos y' 


(2.4) 


Theorem 2.2. The Green function G 77 given by \2.3\ ) has the integral representation 

Jo(\x\rJt 2 — 1) sinh \x\t 


G*{x,y) = S(x)Q(—y) + — sin y 


(cosh xt — cos y) 


-dt 


(2.5) 


where Jq{z) is the Bessel function of zero order. 
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Corollary 2.2. The Green function (2.3) can be represented in terms of the Green function G\ 
as follows: 

G*{x,y)=GZ{x,y)-r [°° y) dt, (2.6) 

J\x\ 




where J\(z) is the Bessel function of first order. 

Remark 2.1. The problem LV = 0 with boundary conditions V(x, 0) = 0, V(x,n) = <p(x) 
can be reduced to the problem (2.1 by the substitution y = ir — y. Using the linearity of L we 
obtain the solution of the problem LV = 0 with the boundary conditions V(x, 0) = tp(x) and 
V (x, 7r) = <p(x) as the sum of solutions of the problems in which one of the boundary conditions 
is a function equal to zero. 

2.2 Green function: construction 


In this subsection we derive the formulas (2.3) for the Green function. To this aim we introduce 
the interaction potential P n (x,y\u) defined as the limit 

P*{x,y;u) = lim P*(x,y;u), 

£—>■0 

The function Pf (x, y, u ) is the solution in the distribution sense (see, e.g. H) of the boundary 
value problem 

LPf = 0, 


P?\ y =o = e-< x ~^e{x - u), 

PI lv=* = 0 - 


We dehne the Green function as the partial derivative in x of P 71 

dP v {x,y\u) 


G(x -u,y) = 


dx 


(2.7) 


( 2 . 8 ) 


Let us consider the Fourier transform in x of the potential Pf: 

v e {v, £, u) = [ e %ix Pff (x, y\u)dx 

Jr 

It solves, as a function of y , the two-point problem 

d 2 "’ 


d y 2 


- (r + r 2 )v e = 0, 


v e (0;£,u) = f^, 
v E (Tr,£,u) = 0 . 


(2.9) 


The solution has the form 


r (V,t,u) = 


e l( - u sinh ^(7r — y)\jr 2 + 


(W 


r 2 + 


(—if + e) sinh 

Making the inverse transform, we obtain the explicit formula for the potential Pf: 

gj£(u-x) /"(7T — y) yjr 2 + £ 2 ) 

P E (x,y,u) = — / --—7— s. 7 d£. 

2n -'K (—+ e) sinh (vr-y/r 2 + £ 2 J 

The integrand is an analytical function in the whole complex plane except zeros of functions 

f + ie = 0, sinh ^7 Ta /^ 2 + r 2 j = 0, 


( 2 . 10 ) 
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Figure 1: Integration loop for the strip case 


that is except the points 


£ = —ie, = ± \Jk 2 + r 2 . 


Integrating along the loop depicted in Figure [I] and applying the Jordan lemma we obtain 
the infinite sum representation of (2.101 


P?(x,y,u) = iQ(x - u)^Res$ e (£+;a;,y) 

k =1 

+iQ(u - x ) -Res$ e (-«e;a;,y) - ^ Res $ s (^;x,y) . 


fc =i 


Here 4> e (£;a;,y) is integrand function in (2.101. Computing the residuals, we get 


iks\n(ky)e^ x u )^ r2 + k2 


sinh ((7T — y)r ) 


Res$ £ (^;a;,y) = =f ‘" kJ “ M "'7' 2 , i 2 \ -> Res x, y) = 

n(r z + k z ) smh(7rr) 

Letting s —> 0 and calculating the partial derivative <9P 7r /dx, we can obtain the result. 


2.3 The general elliptic operator 

Now we extend our formula for the strip of width l and more general operator 


d 2 V n d 2 V 2 

F| y =o = <p{x), l£l, 
V\ y=v = 0, l£l, 


d 2 V 
dy 2 


dV dV 2ta n 

ai— + a 2 ^- r v = 0> 

ox ay 


(x,y) G inttf, 
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where the coefficient p satisfies the condition \p\ < 1. The solution of the problem (2.31 also 
allows the representation V = ip * G l where the Green function G l is given as follows: 


G\x,y) = 6(s{x,y))<d(-y) + 


— R(k)\s(x,y)\ 


(1 - 


e 2 


k =0 v ' 


where 


s(x,y) = —y~x, /3=vv:— 

(J 2 2(1 - P“)cr 1 \ a 2 <T 1 ) 


R(k) = 


1 


2(1 — p 2 )<J\ V a 


a i n P a l a 2 
2 ^ 

0 - 10-2 


+ + 4(1 - p 2 )r 2 + 

a 2 


4(1 — p 1 )u\'K’ 1 k 2 

I 2 ' 


The proof of uniqueness and existence of the solution and the construction of the Green function 
is a straightforward extension of arguments given in Section [4] for the Laplace operator. 


3 Problem on the half-plane 


3.1 Existence and uniqueness theorem 

In this section we study the case of the half-plane II 00 = lx R + . Let us consider the elliptic 
boundary value problem for the operator LV = AV — r 2 V: 

( LV(x,y)= 0, (x, y) £ int II 00 , 

< V\ y =o = <p(x), ieR, (3.1) 

[ V\y = +00 =0, 


Theorem 3.1. Let H x be the space of Holder func tions of order A > 0 and let p be a bounded 
function from H x .Then the solution of the problem (3.1) exists in the classical sense, is unique, 
and allows the representation 


V(x,y)= [ <p(x)G°°(x - u,y)du. 

JR 


(3.2) 


where the Green function G°° has the form 


G°°(x, y) = 6(x)Q(-y) + ± 


2 W o \J £ 2 + r 2 


(3.3) 


Corollary 3.1. The Green function G°° can be represented in terms of the modified Bessel 
function K\ ( z ) of the first order 

G°°(x,y) =S(x)Q(-y) + — /T K i (r\/x 2 + y 2 ) . (3.4) 

Tr^/x z + y z v 7 

Proof. Use the following identity for the Bessel functions (see [3], 3.914) 

- , , A 2 +/3 2 ) , Re/3, Re 7 , a > 0. (3.5) 

\/fl + P v 


irsin(aa;)e PjP+x 2 cry 

dx = — , n 1 


/ 0 y/y 2 + x 2 


□ 
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3.2 Construction of the Green function 


For this domain we use the similar construction of the interaction potential P°°(x,y,u) as the 
limit of solutions P^°(x, y, u) to the boundary value problems 



(3.6) 


Let us consider the Fourier transform 



JR 

Then v e as a function of y solves the boundary value problem for the ODE 



(3.7) 


It can be expressed explicitly: 


e i£u-yy/£ 2 +r 2 



As in the previous case we calculate the inverse Fourier transform and obtain the formula for the 
potential P e °°: 



Integrating along the loop depicted in Figure [2] and applying the Jordan lemma we get the 
following representation for P£°: 



iQ{x — w)Res (<j>(— ie)). 


In contrast to the strip case we have only one residue £ = — ie and two branch points £ = ±ir. 
Letting e —► oo and computing the partial derivative with respect to x, we obtain the needed 
formula for the Green function. 

3.3 The general elliptic operator 

As in the previous section we generalize the results for a more general operator 


f V\ y= o = <p(x), i£t, 
T|y=7T =o, l£l, 


I 



(3.8) 


where p < 1. The solution has the representation V = <p* G°° where G°° is the Green function 



where 
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Figure 2: The integration loop for the half-plane case 


4 Proofs 


4.1 Proof of Theorem 12.11 

We split the arguments into four parts. We prove step by step the following propositions: 


1. The convolution product (2.2) is a continuous and bounded function in domain IP for all 
bounded boundary functions ip and from the space H x . 


2. The convolution product (2.2) is a twice differentiable function in the interior of the domain 
IF for all bounded boundary function ip from H x . 


3. The convolution product \2.2\ is the solution of the problem (2.1) 

4. The solution of the problem (2.1) is unique. 


4.1.1 Proposition 1 


We check the absolute convergence of integrals in the convolution (2.2), implying that the con¬ 
volution is continuous function. Recall the formula from (2.2) 


‘P* G " = - ¥>(«) 


k=0 


ksm(ky)e -\ x ~ u \ Vk2+r2 
\/k 2 + r 2 


du + Q(—y ) / S(x — u)<p(u)du. 


Representing the first summand as the sum of integrals taken over domains (—oo, x] and [x, — oo) 
and using the change of variables x — u = £ in the first integral and u — x = £ in the second, we 
get that 


I nOO 00 

G = - (<p(x + 0 - ip(x - 0) 55 : 

t'o 


(f * (jr = 


ksin(ky)e 


-^y/k 2 -\-r 2 




-d£ + Q(-y)ip(x). 
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The following chain of estimates, where T(i>) is the Euler Gamma function (see pQ) and C denotes 
constants varying from step to step, completes the proof: 




d£ < 


| sin(fcy)| < 1; 

|^(x + C)-¥>(x-C)| < C£ x . 


A ( < 


< 




ke~^ k2+r2 


k—0 


Vk 2 + r 2 


df, < 

k 


f n °° x u ~ 1 e~^ ,x dx = 

JO ’ 

Re/r, Rez' > 0. 


< 


— ^(1 + 0 E , 2 ^,2^1+A/2 — ^ Z-^ fcl + A/2 

k=0 k —0 


^E 


i 


< oo. 


4.1.2 Proposition 2 

We consider the Dirichlet problem for the equation LV = 0 in the disk P centered at the origin 
with radius p. 

LV{x,y) = 0, (i,!/)e!l, 

Plan = v{x,y), 

First we prove the following lemma: 


(4.1) 


Lemma 4.1. The problem (f-1) has a classical solution for any continuous boundary function 
tp{x,y). 


Proof. In the polar coordinates x = p cos 6, y = p sin 9, the problem (4.1) for V(p 1 0) has the 
form 

d 2 V . 1 dV 1 d 2 V _ 


dp 2 p dp p 2 d9 2 

v(p,e)=m- 


(4.2) 


where tp(®) = p{x,y)\(x,y)&an- We separate the variables and solve the Sturm-Liouville problem 
V( P , 0) = f2 Pn(p)Qn(9), + K/P r2Pn ~ Q " 


n—0 


Pulp 2 

It is well known that the solution can be represented as follows 

tf(Odf 


= = n 2 , n> 0. 


v(p, 0 ) = 


it I 0 {rp) J-„ 

- — ^ j" ^sin n8 J sin (nl;) ip (£) df + cos n6 j cos (n^)if(^)df 


where I n (z) is the modified Bessel function (see un¬ 
it is easy to show that I n (z) has the following properties: 


In (z) f Z\ n -0+1 (~) 


In (Z) 


(!) • 


< V 


0(5) 


Oi) r 


In— 1 0) ^ ~ \ n 2 

0.(5) - \z) 


(4.3) 


Indeed, using the definition of the modified Bessel function I n (z ), see jT|, we have: 


T , v'' 00 1 ( z \ [ n+2k v^°° 1 (z\ 

l n (z) _ 2-/k=0 k'.(n+k)\ \2> (Z\ n 2—/k—0 k\(n+k) ! V2/ 


0(5) V°° _ - _ O) 

V 2-/k=0 k\(n+k)\ V 2 / 


n+2fc — ' ' 


(!) 


E oo 

k =i 


2k 

,2 k< I' 


k=0 k!(n+k)! \2/ 


O' 


8 































l n+l\ 


Z) EZ 


0 fc!(n+fc+1)! 


(!) 


n+l+2fc 


E OO 1 ( z \ 

k= 0 fc!(ra+fc)! V 2 / 


2k 


4,(5) v°° _i_fiy 

2^k=0 k\(n+k)\ \2> 

In— l(© k\(n+k — 1)! (2) 


< - t- 


n— 1+2/c 


( 4 ) 

Vz/ 2 v°° 1 tit 

Z-^k—0 k\(n-\-k)\ \2 ) 

E oo 1 / 2: \ 

fc=0 fc!(n+fe-l)! UJ 


© 


< - ", 


2fc 


, n+2/c — 


© 


< - - 


4,(5) V°° _i_(© 

2-Jk=0 k\(n+k)\ V 2 / 

From here we immediately obtain the bounds for V ( p , 9) and 

d 2 V(p,9) 


E OO 1 / z \ 

k—0 k\(n+k-l)\ V 2 / 

d 2 v( P ,e ) 


2k — 


© 


< - 


co / \ n 

\v{ Pl e)\ < ( 4 J < Const, 

n=0 


90 2 


<^E 


n—0 


56»2 


n 2 ( 4 ) < Const. 

P; 


For the derivative dV/dp we use the properties of the modified Bessel functions, see IT]: 

d Tl 

I n -i(rp) + I n+ i(rp) = 2—— ( I n {rp )), I n -i{rp) - I n+ i(rp) = 2 — I n {rp). 

rap rp 


Thus, the derivative dV/dp is bounded because 

+ £ 


dV 


In-i(rp) y- In+1 {rp) 

fb 7 «(rp) / n ( r © 


< Const. 


These relations mean that V, dV/dO , and dV/dp are continuous functions. The arguments for 
the second derivative d 2 V/dp 2 are similar. □ 

Now we can complete the proof of Proposition 2. We take an arbitrary point ( x,y ) from 
intn 71 ’ and consider a disk centered in (x, y) and contained in intll 71 '. Due to Proposition 1 the 


convolution (2.2) is continuous on the boundary of this disk. By Lemma 4.1 (2.2) it is twice 


differentiable at any point of its interior. Hence, the convolution (2.2) is a twice differentiable 
function in the interior of n ir . 


4.1.3 Proposition 3 

First, we check the boundary conditions 


V(x,y)\ y =o= / (p(u)G 7r (x - u, 0)du = <p(x), V(x,y) \ y =n = / (p(u)GF(x - u, n)du = 0. 

JR JR 

It is easy to show that the Green function is a weak solution, i.e. we understand the equality 
LG lr = 0 in the following sense: 

LG" = 0 ((G© Lz) n = f f G n Lzdxdy = 0, Vz £ G 2 , VH C W \ , 


where C 2 is the class of twice continuously differentiable finite functions (the class of test func¬ 
tions) and H is a compact sub-domain. Hence, for the convolution ip * G 7r we have 

(<£> * G© Lz) n = (( (<p * G© Lzdxdy = <p * ff G n Lzdxdy = 0. 

J jo J Jo, 

It is well known that the elliptic operator is self-adjoint. Therefore 

((L{<p * G©, <p) n = (<p* G© L*z) n = © * G© L*z) 0 = (ip* G*,Lz) n = 0. 

Due to Propositions 1 and 2 function the L(ip*G n ) is continuous and, therefore, bounded on the 
compact D. Hence, L(tp * G© = 0 because fi is arbitrary. 
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4.1.4 Proposition 4 


Suppose that we have two different bounded functions V and V, both solving the problem (2.1). 
Their difference V — V is a bounded function solving the homogeneous problem 


LV 0 (x,y) = 0, (x, y) e int n 71 

Vo(x, 0) = 0, x € R, 

Vq{x, 7r) = 0, ieK. 


(4.4) 


Lemma 4.1. The solution of the problem (4-4) is the infinite sum 


OO 

V 0 (x, y) = J2 [A k e V ^+ Wx sin ky + B k e~ V ^+ Wx sin 


k =0 


(4.5) 


It is easy to show that all summands in (4.51 are unbounded functions except zero (the 
eigenfunction corresponding to k = 0). Hence, V = V. 


Proof. We separate the variables Vo(x,y) = X{x)Y{y) and solve the Sturm-Liouville problem: 

Y" _ Y —V" 

YX" + Y"X — r 2 XY = 0, -—-= —— = A 2 . 

X Y 

Spectrum is A = X k = k, k £ N and the eigenfunctions are X k (x)Y k (y) = e ±vV 2 +fc 2 :c s jn 


4.2 Proof of Theorem 12.21 

We put £ = x — u and consider the part of the Green function G 7T given by the series: 


R(€,y,r) 


oo 


E 


fcsinfcy e? ^2+y2 
\fk 2 + r 2 


(4.6) 


In the case r = 0, see [3] 1.445.1, we have: 

siny dR{f,y, 0) 


R(Z,v, o) = 


2 (cosh£ — cosy) ’ 


df 


= T(£,y, 0) = - 


sin y sinh f 
2 (cosh£ — cos y)“ 


(4.7) 


Using the formula 


see 0, 6.646(1)), we have 


/ OO 

e~ k ^J 0 (fr\/1 2 — l)dt = 


e ~iVk 2 +r 2 

fVk 2 + r 2 ’ 


°° /-oo _ 

!?(£, y,r) = £ ^ /c sin fcy / e~ k ^Jo [fr\/1 2 — lj dt. 
fc=o 


The change of the integration and the summation yields 

R{€,y,r) = Jo (^ry/t^Yj T(£t,y, 0)dt = -J J 0 (^rVt 2 - lj d (R(£t, y, 0)). 


Using the second formula in (4.7), we finish the proof. 
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4.3 Proof of Theorem 13.11 

The proof consists of four parts: 


1. The convolution product (3.2) is a continuous and bounded function in the domain II 00 
for all bounded boundary functions <p from H x . 

2. The convolution product (3.2) has the continuous second-order partial derivatives in 
IT^/SII 00 for all bounded boundary functions <p from H x . 

3. The convolution product ( |3.2|) is the solution of the problem (3.1) 

4. The solution of the problem (3.1) is unique. 

We omit the proofs of Propositions 2 and 3 as the arguments are similar to those in the 
previous case. 

4.3.1 Proposition 1 

We make changes in the variables similar to those used above and get that 


G°° * p = - 


sinl/b/lp 

(ip(x + 77) - <p(x - 7?)) -- - drjdt; + G(-y)ip(x). 

V?’ 2 + f 2 


Using the bounds relations we infer the absolute convergence of the integral: 

\fr' 1 + £ 2 


< 


J Js.1 \/r 2 + £ 2 


< 


|sin(fey)| < 1; 

+ 0 -kp(x-C>\ < ( - 


< 


/ 0 oo ^- 1 e-^x = ^, 
Re fj,, Re v > 0 


< CT(1 +A) 


n r-v —- 

l V r +£ 

ddt 


drjdti < 


'0 (r 2 + £ 2 ) 1+V2 


< C 


dl; 


^l + A/2 


< c. 


4.3.2 Proposition 4 


As in the prev ious case we suppose that there are two different bounded solutions V and V of the 
problem (3.1). Then we have the homogeneous boundary problem for the difference Vq = V — V. 


A Vo - r 2 Vo = 0, (x, y) G int II 00 , 

Vo(a;, 0) = 0, i£l, 

Vq(x, 00 ) =0, 


(4.8) 


We show that the solution of (4.8) can be represented in terms of the modified Bessel functions. 
After that we use the asymptotic of modified Bessel functions and show that only the zero 
function solves the problem (4.8) in the class of bounded functions. The proof of the lemma 
below completes the proof of Proposition 4. 


Lemma 4.2. The solution of the problem (4-8) is the infinite sum 

OO 

V 0 (x,y) = ^ ~2(A n I n (ru ) sin ny + B n K n (rv) sin ny). 


(4.9) 


n =0 


This sum does not have bounded summands for n > 1. 


Proof. In the polar coordinates x = v cos p and y = v sin ip we have the following Sturm-Liouville 
problem in separated variables V = N(i/)$>((p): 

N" + N'/v - r 2 N 4" 

N/v 2 ' 
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For the phase component $ we have the spectral problem 


$" + A 2 4> = 0, 

$(0) = 0, (4.10) 

$(7r) = 0. 


The eigenfunctions are <bfc(</?) = sinfc^> and A = A k = k with k € N. For the radial component 
we have the modified Bessel equation 


M' 1,2 

K + — - r 2 N k - ~^N k = 0 . 

v 


Therefore, the solution of the problem (4.8) can be represented in the form from (4.9). The 
modified Bessel functions have the well-known asymptotic, see [T|, 

K n (z) ~ oo, l n (z)~0, (z-^-0), K n (z) ~ 0, I n (z) ~ oo, (z-too). 

Hence, in the sum we do not have bounded summands for n > 1. □ 


Acknowledgement. The research is funded by the grant of the Government of Russian Feder¬ 
ation n°14.12.31.0007. 
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